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Abstract. By combining the first year data from the Supernova Legacy Survey (SN 
LS) and the recent detection of the baryon acoustic peak in the Sloan Digital Sky 
Survey data, we are able to place strong constraints on models where the cosmic 
acceleration is due to the leakage of gravity from the brane into the bulk on large 
scales. In particular, we are able to show that the DGP model is not compatible with 
a spatially flat universe. We generalize our analysis to phenomenological toy models 
where the curvature of the brane enters into the Friedmann equations in different ways. 
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1. Introduction 

Observations of Type la supernovae have led us to understand that the universe started 
to accelerate at a redshift of ^ ~ 1 [TJ, |2j [HI - This suggests the presence of some energy 
component which redshifts slowly and which had the same density as matter relatively 
recently, z < 1. The desire to explain this behaviour has led to the development of 
a number of models such as quintessence where the energy density corresponds to the 
potential energy of some scalar field. 

At the same time, models with large extra dimensions into which only gravity 
can propagate have become popular which suggests that there may be short distance 
corrections to gravity j2J El- When the Friedmann equations are solved in these 
higher dimensional backgrounds, the expression for the Hubble parameter changes to 
if 2 = G(p + kp 2 ) where k is related to the tension of the brane upon which the standard 
model particles reside [Zj. However, modifications of this form do not give rise to the 
kind of behaviour required to explain the observed acceleration. 

A related theory is that due to Dvali, Gabadadze and Porrati (DGP), which presents 
a class of brane theories where the theory remains four dimensional at short distances 
but crosses over to a higher dimensional theory for gravity at some very large distance 
(Hj. In this model, gravity leaks out into the bulk only at large distances and the resulting 
modifications to the Friedmann equations can give rise to universes with acceleration 

[gnu. 

In section 12] we will briefly present the DGP model and then we will show in section 
El that is it marginally disfavoured by observations of supernovae and baryon oscillations, 
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at least in the case of a spatially flat universe. In section 0] we will try to generalise 
the 1st (tt) Friedmann equation obtained in the DGP model to see if the data indicates 
what kinds of modifications to the Friedmann equations might explain the data. 

2. The DGP model 

The DGP theory starts with the idea that standard model gauge fields are confined 
to live on a (3+l)D brane residing in an infinite (4+l)D bulk, with different scales of 
gravity on the brane and in the bulk [Hj. The gravitational part of the action is therefore 
given by 

S =^J d 5 zV=^R 5 + J^J d 4 x^g~ 4 R A (1) 

where p 2 = Mpf and k 2 = M^ 3 . The Friedmann equation which comes from this 
gravitational action coupled to matter on the brane take the form JT] 



rv r\j I T o \ iX 
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where e = ±1. For our purposes it is important to note at this point that the corrections 
to the normal equations involve the extrinsic curvature of the brane itself, in other words 
H 2 + ka~ 2 . Later when we try and generalise this equation to compare with the data, 
this observation will loosely motivate the form of our generalisation. If we set e = 1 we 
can re-arrange things so that the Friedmann equation takes the form 

2 
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where the length scale r c is defined to be 

± = tL 

2r c k 2 

and r c sets a length beyond which gravity starts to leak out into the bulk, modifying 
the Friedmann equations. The relationship between the Hubble expansion parameter 
today Ho and that at a redshift z is then given by 

H(z) 2 



( 4 ) 



H 2 



H'{zf = n k (i + zf + [y/n M {\ + zf + n rc + v^CT ( 5 ) 



where we use the usual notation Q x = p x /p cr it where p x refers to the energy density in 
component x today and p crit = 3Hq/8ttG. The new unusual parameter Q rc takes the 
form 

,4 



/' 

'0 

and the constraint equation between the various components of energy density at z = 
is given by 

Q M + n k + 2 V ^T V / 1-^ = 1 (7) 
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So that, in particular for a flat universe with k = we have the relation 



It has been shown that by setting the length scale r c close to the horizon size, this extra 
contribution to the Friedmann equation leads to acceleration which can in principle 
explain the supernova data [01 E3] . In the next section we will see if this theory can fit 
with the latest observations. 



3. Confronting the model with cosmological data 

Measurements of the brightness of Type la supernovae as a function of redshift are 
sensitive to the cosmological model via the integration over expansion history in the 
expression for the luminosity distance (c = 1) 



H ^\Ch\ V "Jo H'{z) 

where the function S(x) is defined as sin(x) for < (closed Universe), sinh(x) for 
> (open Universe) and S(x) = x, and the factor is removed for the flat 

Universe 

In the DGP model, H(z) is then given by Eq.(jSJ. Very recently, the SNLS 
collaboration published the first ear data of its planned five year survey (Astier et 
al, 2005) ^2j. The data set includes 71 high redshift Type la supernovae (SNIa) in 
the redshift range z=[0.2,l] and 44 low redshift SNIa compiled from the literature but 
analysed in the same manner as the high-z sample. Thanks to the multi-band, rolling 
search technique and careful calibration, this data-set is arguably the best high-z SNIa 
compilation to date, indicated by the very tight scatter around the best fit in the Hubble 
diagram and the careful estimate of systematic uncertainties in jT2j. In our analysis (as 
in [T2j) we fit for cosmological parameters with the prior from the baryon oscillation 
peak detected in the SDSS Luminous Red Galaxy survey (LRG) of Eisenstein el al 
(2005) [in!, g iven b Y 



M 



1 / n r C Z1 dz 

,s MoT 
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0.469 ±0.017, (10) 



where E(z) = H(z)/H and Z\ = 0.35. Furthermore, S and are defined as in Eq.(JUJ). 
The quoted uncertainty corresponds to one standard deviation, where a Gaussian 
probability distribution has been assumed. Figure El shows the allowed regions in the 
Q rc -flM plane in the DGP model from each of the two data sets and the combined 
limits. Also shown in the graph is the curve corresponding to flatness, i.e. Qk = 0. 

The allowed region of the VL Tc -VLm plane deduced from baryon oscillations forms a 
band whereas the allowed region from supernovae luminosity measurements forms an 
ellipsoid. The fact that the baryon oscillations constrain models to lie on a band reflects 
the fact that they have only so far been detected at two redshifts, z = 0.35 in the SDSS 
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data and z ~ 1100 in the CMB data. This gives us one ratio vs. the many ratios 
which come from the many type la supernovae. However, the power of this approach 
lies in the fact that the band of allowed parameters from baryon oscillations is almost 
perpendicular to the major axis of the ellipsoid coming from the supernovae data. This 
occurs in a very similar way in the conventional Q^^m plane. The use of the baryon 
oscillation prior therefore allows us to constrain model universes much more effectively. 

We conclude that the fit of the DGP model to the supernova data along with a 
prior from the baryon oscillation results is incompatible with a flat universe. Given the 
evidence for the latter from the CMB experiments, this model is clearly disfavoured by 
the current observational data. For comparison, in FigureOlwe show the fit to the same 
data to a dark energy model generically described by an equation of state parameter, 
w(z) — Wq + Wi ■ z, where a flat universe is assumed. The A model (wq = —l,w\ = 0) 
gives an excellent fit to the data. All the fits were performed with a modified version of 
the cosmology fitter (SNALYS) in the SNOC package [Hj. 



4. Generalisations of the DGP model 



We would like to consider generalisations of the DGP model in order to see if the data 
favours any particular modification of the Friedmann equation. (This is in the spirit of 
15 j and 16J.) The extra term in the Friedmann equation (J2J) comes from the derivative 
of the scale factor in the 5th dimension, which is set by the Israel branching conditions 
across the brane and therefore by the extrinsic curvature on the brane. We will therefore 
assume the extra term enters in the following form 

fl,+ ?r=$( B,+ ?)-T«- < n) 

where now a is a free parameter and n 2 is now related to some phenomenological mass 
scale M 

nl = M a ~\ (12) 
The expression for the Hubble expansion as a function of redshift is then 

or equivalently 

H( 7 \ 2 , //7(V1 2 \ a/2 

-A±L-n k (i + z) 2 -2^i c - n k (i + z) 2 ) = n M (i + z) 3 (u) 

where now we have generalised the expression for the DGP brane world dark energy 
density fraction 

Jn7 '= (15) 

which means that the constraint between the various densities at z = is given by 

Q M + Q k + 2v^T(l - VL k ) a ' 2 = 1 (16) 
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Figure 1. The solid curves show the allowed parameter regions in the S1m — plane 
at the 68, 90 and 99% confidence level. Solid thin contours correspond to the the first 
year SNLS data in Astier et al (2005) ^3 ; the- dashed lines show the corresponding 
regions from the baryon oscillation peak in Eisenstein et al (2005) JH]- The coloured 
contours indicate the result of the combination of both data-sets. The thick solid lines 
indicate the expected relation between Qm and f2 rc in a flat universe. 
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Figure 2. Allowed 68, 90 and 95 % confidence regions in the wq,wi parameter 
space for the first year SNLS data with baryon oscillation prior in a flat universe. 
The linear expansion w(z) = wq + w\ ■ z has been assumed. The best fit solution 
(uif) = —1.07, wi = 0.14) is also shown. 



In order to make sure there are no pathologies here, we need to ensure that 

n k <i-n M ■ (n M >o) (17) 

then if we write 



'H(z) 2 ^ k 



m ami 



we can then solve the following equation as a function of redshift 

h 2 -2y^ c h a = tt M (l + z) 3 (19) 



and use the solution to obtain the luminosity distance. Since a is no longer a = 1 we 
now obtain the solutions to this equation numerically. 

Figure E| shows the allowed values in the o>-VLm parameter space if a flat universe 
prior is imposed. The best fits are centered close to a = 0, which corresponds to 
a constant in Eq.((T3J), i.e. equivalent to A. The full data set of the SNLS, ~700 
spectroscopically confirmed Type la SNe, will soon provide further constraints in the 
possible values of the exponent a. 
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Figure 3. Allowed parameter regions in the (x-VIm plane at the 1 a (one parameter), 
68, 90 and 95% confidence level for a flat universe. 

5. Conclusions 

The most recent high-redshift supernova data has been used to test the DGP brane 
model of cosmology where dark energy is replaced by an extension of the Friedmann 
equation to account for the leakage of gravity to a large extra dimension. Using a prior 
from the SDSS baryon oscillation results we find that the this model is incompatible 
with the data for a flat universe. Generalizations of the model including an extra 
term H a /r^~ a yields —0.8 < a < 0.3 (1 a) i.e consistent with a = 0, the situation 
corresponding to a non-redshifting cosmo logical constant. 

This is not completely surprising if we look at the generalised Friedmann equation 
(jlljl since during the matter dominated era, it is straightforward to show that the 
equation of state w = — 1 + f jT^j. Observational constraints on w already push us close 
to w ~ — 1 (see Figure EI) so we see that we require some small value of a in order to 
match the observations at this time. 

Finally, we note that it may be interesting in the future to see if it is possible to 
obtain a good fit to the data in models which give rise to different modifications to the 
Friedmann equations. 
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